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Abstract 

Three variants of mean field methods for atomic and nuclear reactions are compared 
with respect to both conception and applicability: The time-dependent Hartree-Fock 
method solves the equation of motion for a Hermitian density operator as initial value 
problem, with the colliding fragments in a continuum state of relative motion. With 
no specification of the final state, the method is restricted to inclusive reactions. The 
time-dependent mean field method, as developed by Kerman, Levit and Negele as well as 
by Reinhardt, calculates the density for specific transitions and thus applies to exclusive 
reactions. It uses the Hubbard-Stratonovich transformation to express the full time- 
development operator with two-body interactions as functional integral over one-body 
densities. In stationary phase approximation and with Slater determinants as initial and 
final states, it defines non-Hermitian, time-dependent mean field equations to be solved 
self-consistently as boundary value problem in time. The time-independent mean field 
method of Giraud and Nagarajan is based on a Schwinger-type variational principle for 
the resolvent. It leads to a set of inhomogeneous, non-Hermitian equations of Hartree- 
Fock type to be solved for given total energy. All information about initial and final 
channels is contained in the inhomogeneities, hence the method is designed for exclusive 
reactions. A direct link is established between the time-dependent and time-independent 
versions. Their relation is non-trivial due to the non-linear nature of mean field methods. 



1 Similarities and Differences 



After the success which mean field methods have had for bound state problems in various 
fields of physics, it was only natural to try the mean field concept for scattering states as well. 
The original attempt in this direction is the time-dependent Hartree-Fock method (TDHF) 
developed about 20 years ago [1]. In this method one solves the equation of motion for the 
one-body density operator p = pit), 

d 

ih — p= [h, p] , (1.1) 

with an initial condition for p, 

P (U) = pi , 

describing an unbound state of relative motion of the colliding fragments 
p can be represented in just one basis of single-particle states ip m (£), 

N 

P= E \^m{t) >< i> m (t) | (1.3) 
m = 1 

for N particles and with %p rn (t) orthonormalized. The single-particle Hamiltonian h is Hermi- 
tian and has standard Hartree-Fock structure, h = t + u with mean field potential 

N 

U= E < -^m(f) \v\ • 1p m (t) > , (1.4) 
m = 1 

where the dots represent wave functions of some arbitrary single-particle basis. The single- 
particle states ip m (t) are determined from the time-dependent Hartree-Fock equations 

d 

in—\^ m {t) > = h\il> m (t) > (1.5) 

with initial conditions 

ipmiti) = Xm for m = 1, 2, ... N . (1.6) 

Standard manipulation of eq. (1.5) and its Hermitian adjoint leads to eq. (1.1), with initial 
condition (1.2) corresponding to (1.6). In practice, the non-linear equations (1.1) and (1.5), 
respectively, are solved by iteration: A set of iV occupied orbitals is used to compute the poten- 
tial u. With the corresponding solutions of (1.1) or (1.5) one recalculates u till self-consistency 
is reached. As each iteration step requires as input only the set of iV functions ip m at fixed time 
t, the TDHF method as initial value problem is "local" in time. 

From the densities Pbif, t) according to (1.1) and (1.2), taken at various classical impact pa- 
rameters b, one may calculate a classical cross section. Although it can generate nice snapshots 
of the density distribution during the scattering process, the method has two problems: First, 
as initial value problem it can at best handle inclusive reactions, since no specification of a 
final channel enters the formalism. The second, more serious problem concerns "spurious cross 
channel correlations" : Starting from the determinant \& (t«) corresponding to pi the method 
generates a single determinant ^ (t). Due to the non-linear nature of the TDHF equations, 



(1.2) 

. The density operator 
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this determinant \I/ (t) continues to vary as time goes to infinity. Hence, if this wave function 
is expanded in an orthogonal set of channel wave functions, the expansion coefficients will not 
be constant asymptotically. Thus an S-matrix, constructed by projecting the TDHF wave 
functions onto channel wave functions, would not be constant in time [2]. 

The time-dependent mean field method (TDMF) [3, 4] also uses an equation of motion like 
(1-1), 

d 

ih — p= [h, p] , (1.7) 
however, there are two important differences: 

1. The density operator p of (1.7) is expanded in two sets of mutually biorthogonal single- 
particle functions, 

p = j2 1 (*) >< $™ (*) I ( 18 ) 

m=l < 4>m(t) \ lpm(t) > 

and the Hamiltonian h = t + u with mean field potential 

m = 1 

is non-Hermit ian in general. 

2. Eq.(1.7) has to be solved self-consistently as boundary problem in time t, fixing 

p(U) = pi and p(t f ) = p f . (1.10) 
The single-particle functions ip m (t) are obtained from 

ih^- t \t/j m (t) > = h\jjj m (t) > (1.11) 
by forward propagation of initial wave functions 

^mik) = Xm for m = 1, 2, ... N . (1.12) 
Analogously ■0 m (t) results from 

- ih— <4> m (t)\ =<4> m (t)\h (1.11') 

by backward propagation of 

^m(tf) = X m for m — 1,2, ... N . (1.12') 

Combining (1.11) and (1.11') in the usual way, one obtains (1.7) with conditions (1.10) from 
(1.12), (1.12'). One also proves easily that 

^ <i) m \i) n > = , (1.13) 
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hence if ip m , ip n are chosen biorthogonal at t — ti, they will remain so at any time. In general, 
ip m and ifj m will be not complex-conjugate to each other, and p will be non-Hermit ian. Only 
when x m is generated from Xm by the mean field time development operator Uh(tf — ti), we 
have ip m (t) = ip m (t) in which case TDMF reduces to TDHF with just one set of single-particle 
functions. Equations (1.7), (1.10) as well as the coupled equations (1.11), (1-H') together with 
(1.12), (1.12') constitute a boundary condition problem in which the time plays a role similar 
to that of the spatial coordinates [5] . 

In practice, one may try to solve this problem self-consistently in analogy to the static Hartree- 
Fock problem by iteration [5] which involves an initial guess of the single-particle functions 
ip m (t), ip m (t) for all times t between ti and tf. In this sense, the TDMF boundary condition 
problem is highly "non-local" in the variable t. Since initial and final states are taken care of, 
the method is able to describe also exclusive reactions: Each quantum process, leading from 
given initial state x to some final state x > nas its own time-dependent mean field assigned. It 
has been proven [2] that, within the framework of TDMF, an S-matrix can be defined which 
becomes asymptotically constant. The problem with the TDMF approach lies in the above 
non-locality in time, for which there seems to exist no practicable algorithm for actual numer- 
ical calculations of (3 + 1) -dimensional systems. 

Some comments on the lack of Hermiticity of p and of h (p) are in order. This non-Hermitian 
structure is of a special type which preserves the usual properties of a single-particle density 
matrix connected to a Slater determinant. From the explicit expression (1.8) we have immedi- 
ately 

Trp = N and p 2 = p , (1.14) 

hence particle number is fixed and p is a projector on the space of occupied orbitals. Moreover, 
| ip m > and < ip m \ are right and left eigenstates of p, 

p | ip m (t) > = N m | ip m (t) > ; < ^ m (t) | p = < ^ m (t) | N m (1.15) 

with eigenvalues 




1 for m occupied 
otherwise 



The eigenvalues of h (p) will be complex in general, hence ip m and ip m correspond to quasi- 
particles of finite lifetime. They represent intermediate states which the system passes during 
the reaction process. The S-matrix of the mean field approach will always violate unitarity, 
irrespective of the fact whether h (p) is Hermitian or not: The iS-operator depends on p which 
in turn depends on the initial and final states of a specific reaction as mentioned above. Each 
quantum process x ~^ X nas its own mean field, calculated without reference to any other 
possible process. 

An alternative to TDMF is the time-independent mean field method (TIMF) based on a 
Schwinger-type variational principle [6, 7] as described below in section 2. Like TDMF it 
uses two sets of variational functions and leads to a set of inhomogeneous equations for the 
single-particle functions <p i5 ipi which may be chosen biorthogonal. The corresponding density 
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operator 



N 

a < r / I W > 



(1-17) 



obeys an inhomogeneous equation, 

= 570v -P X ) 5 ?7 = £ - ^-Ffrf^ » ( L18 ) 





if 


0> 


<0' 


0> 



where E is the total energy of the system, 0, </>' are Slater determinants built from the single- 
particle functions ifi, <f i and p x >, p x are "mixed" densities containing all information about 
initial and final states x, X ■ Hence the method can describe exclusive reactions. As in TDMF 
h is non-Hermitian, 



A' 



< -<p'i\ 


V 


1 • <p* > 


< <p'i 1 





h = E z i , (i.i9) 



but now of course time-independent as are the single-particle functions tpi, </^. Their self- 
energies r]i are complex in general due to the lack of Hermiticity of h. As regards its practical 
applicability [8], TIMF is much simpler than TDMF: One has to solve inhomogeneous, complex 
equations of Hartree-Fock type at some given energy E. In this respect TIMF is comparable 
to TDHF: While TDHF is "local" in time, TIMF is "local" in energy. 



In view of the above similarities and differences of time-dependent and time-independent mean 
field methods (see table 1), it appears desirable to establish a direct link between TDMF and 
TIMF. Their relation is non-trivial due to the non-linear nature of mean field methods. 



The paper is organized as follows: In section 2 we shall derive matrix elements of the resolvent 
operator as stationary values of a Schwinger-type functional in mean field approximation. Its 
stationarity conditions are given by the TIMF equations (1.17) and (1.18). To calculate matrix 
elements of the exact time-development operator we use a path integral representation based on 
the Hubbard-Stratonovich transformation (section 3). Ambiguities of the respective auxiliary 
fields a are briefly discussed in section 4, following [9]. The mean field approximation of the 
time-development operator is then obtained by applying the stationary phase approximation 
to the path integrals over the auxiliary fields a in section 5. The result is the boundary 
condition problem for p, eqs. (1.7) to (1.10), of the TDMF method. The link between the 
two methods is obtained in two steps: First, in static approximation of TDMF (section 6) we 
can introduce single-particle propagators g m = (w m — /i(ct)) _1 by Fourier transformation of 
the corresponding time-development operator exp {— | h(a) T}. Second, the Fourier integrals 
over w m are solved in stationary phase approximation on the same footing (section 7) as the 
path integrals over the static auxiliary fields a. Thus the integrals over real variables w m are 
replaced by the integrand at some complex values io° m which are identified in section 9 as the 
self-energies i] m of TIMF after the T-integration has been carried out (section 8). At the same 
time one reproduces the TIMF value of the resolvent. In conclusion (section 9): TIMF turns 
out to be a static approximation of TDMF in semi-classical approximation. 
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2 Resolvent operator in mean field approximation 

The TIMF method solves the stationarity equations of a functional like [6, 7] 

F(tf', *) =< x'l* > + < > - < - > (2.1) 

which read 

(E - H)\*>= \ X >; <*'\(E - H) =< X '\ (2.2) 

with Xi X as initial and final states. The resolvent matrix element between x an d X is then 
given by 

<X\(E - H)- 1 \x> = <x'\y > = <^'\X> ■ (2.3) 
In mean field approximation one takes x, X as Slater determinants, 

N N 

x = A n ; x' = ^ II x'i, (2.4) 

i=l i=l 

and the variational functions VP' are correspondingly restricted to Slater determinants as 
well, 

N N 



</> = AYltpi ; 4 = A I] Pi ■ (2-5) 

i=i i=i 

The antisymmetrizer ^4 is defined as A = (AM)" J2 (~ ) P P- One then has to vary 



p 



F(0',0) =<x'|0> + <0'|x> - <<f>'\(E - H)\<j>> . (2.6) 

To calculate F (</>', 0) and its functional derivatives, one can use the invariance of the above 
four determinants under unitary transformation of their orbitals in order to diagonalize the 
four iV-dimensional matrices < x% I Pj >> < Pi I Xj >> < Pi I Pj > an d < Pi I h \ ipj >. 
Here h = t + u is the single-particle Hamiltonian with -u the Hartree-Fock potential in the 
biorthogonal basis p^Pj- In this representation the stationarity equations of (2.6) read: 

(rji - h)\<pi>= K\Xi> ; < Pi I (Vi ~ h) = < Xi I A- (2.7) 
with complex self-energies 



m = E- " y '/ T + ^Hf^ • (2.8) 



<0' 


H 


0> 


<0' 







h 


Pi > 


<<Pi 


Pi > 



The stationary value of F ((f)' , 0) is then 

^ = < x' I > = < 0' I X > = II A * < Xi I 9i I Xi > = II A i < Xi I 9i I Xi > (2.9) 

i i 

with single-particle propagators 

9i = (Vi- h)- 1 . (2.10) 
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The explicit form of Aj, A^ - which we shall need in the following - is 

x < <t> 1 X > < v'i\Vi> x ' < X I > < ^ I > , 9 1 1 x 

4 <0'|0> ' <<Pi\Xi> ' * <0'|0> ' <Xi\<Pi> ' 1 J 



They are related to each other and to rj of eq. (1.18): From (2.11) and (2.9) one finds 

< 0' I X > \ N1 = ( < X I > 

< 0' I > / ' V < 0' I > 



^ = =nt (2.12) 



Furthermore 



^ 1 = n ^ = n ^ , (2.13) 



as is easily seen: Scalar multiplication of (2.7) by < </?^ | and | y?j >, respectively, gives together 
with (1.8) and (2.8): 

E _ «a\H\<i» <rijfehgi> = Ai <y;ixi> = A ; <x;i^> (214) 

< I > ' < I fi > <¥i\¥i> <¥i\¥i> 

and with (2.11) 

_ < 0' | X > _ < X I > f2 15 , 

^ ~ <0'|0> ~ <0'|0>' ( 5) 

Combining (2.12) and (2.15) then results in (2.13). 

It is now easy to obtain an equation for the density operator 

= » ^><rij , (2 . 16) 

using biorthogonality but leaving the normalization of the single-particle functions open. The 
dyadic product of (2.7) with < ip i | and | ipi > from right and left, respectively, gives 

(rji - h) | ^ >< ip\ | = Xi | Xi >< <Pi I ; I <P* >< <Pi I (Vi ~ h) = | >< xl I ^ . (2.17) 
Dividing by the norm < y?- | >, summing over i and forming the difference finally results in: 

[h, P ] = V(P X ' ~ Px) , (2-18) 

if one defines the "mixed" densities 

p, = £ lll^Al . px = £ I * >< 4 I (2.19) 

and uses equation (2.14). This constitutes the proof of equation (1.18). 
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3 Path integral representation of the time— development 
operator 

If a state x is prepared at time t iy then it will develop with time upto t = tf as 

\X (*/)>= -*i)|x> (3-1) 

with 

[/ (t f - U) = exp |- % - (tf - U) if} . (3.2) 
The probability of finding some state | x' > at £ = is then obtained from the amplitude 

<x'\x(tf) > = <x'\U(t f - U)\ X > ■ (3.3) 

An exact representation of the time-development operator of the many-body problem in terms 
of a time-dependent single-particle Hamiltonian can be achieved by the Hubbard-Stratonovich 
transformation. To this end one rewrites the full Hamiltonian if in terms of operators [3] 



(3.4) 



in arbitrary orthonormal basis as 



(3.5) 

a, 7 1 a,P,f,8 

with matrix elements 

vaw = J J d 3 fd 3 f' V l(r)^(f')v(f,f')^(r) V5 (r) = (a (3\v\l , (3-6) 

t ai = J d 3 f(p* a (f)t(p 1 (f) and K ai = t aj - ^ v *i 



in that basis. Whenever possible, we shall use the matrix notation 



H = K ■ p+ X -p ■ v ■ p (3.7) 

to avoid cumbersome labels and summations. The unphysical self-interaction term in (3.6), 
combined with the kinetic energy t ai into K ai , arises from anticommuting creation and an- 
nihilation operators in the standard form of if such that the density operator representation 
(3.5) follows. We shall come back to this point in section 4. One then employs the Gaussian 
trick in its complex version, 

exp • A ■ &} = VdetA- 1 J -J=j exp j + ^rr • A' 1 ■ x - ib ■ a:} , (3.8) 

to the time development operator (3.2), identifying the density operator p with b. Relation 
(3.8) is valid if matrix A is real, symmetric and invertible. 
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The Gaussian trick (3.8) cannot be applied immediately to the two-body part of H in equation 
(3.5), since the one-body and two-body parts of H do not commute. This non-commutativity 
is treated the standard way by dividing the time-interval (tf — U) in (3.2) into M slices of equal 
length e — (tf — t^/M. Then we can factorize for e — > 0, 



exp {~^h} = e W [- l -eK.p)^{-^_ep.v.p) , 
and we may now apply (3.8) to linearize the second factor of (3.9) in p: 



(3.9) 



exp ' ~2h ep ' v ' 9 



det i^v- 1 ) 



d (y a j 



a — —ea-p 

n 



For each operator p ai we have to introduce a real variable a aj . Choosing 



a — a ■ v 



(3.10) 
(3.11) 



we may also write 



as useful alternative to (3.10). While a has the quality of a density, o has that of a potential. 
Matrix t> of equation (3.6) is symmetric in particle coordinates r, r and thus in the label pairs 
(a, 7) and (/3, 5). Hence v fulfills the symmetry requirement on A in formula (3.8). 

Repeating the above step for each time-interval, labelled by index k in the following, we can 
write for the time-evolution operator (3.2) with (3.10): 

f i M 1 M r % ^ 

x ex P I ^ eY,v(k)-v- 1 -a(k)jl[expt [ --e(K + a(k))-p ] f . (3.13) 

To prepare the next step we note that from (3.10) 

( det (l )-W(n^)exp{^. „-..}. (3,4) 

We may now take the limit e — > 0, M — > 00 such that eM = (tf — U) remains finite, replacing 
in (3.13) 

*/ 

— ► J dt ■■■ . (3.15) 



M 

e E 

fc=i 



Then (3.13) reads as functional integral, using the Trotter formula, 



7 \j2ii% , 



(3.16) 



r . */ I r . tt 

X 6XP I 2~ft / d * °" ^ ' ' °" ^ ' T ' 6XP I ^ / dt ^ K + a (*)) ' ^ 



where T denotes time-ordering and the norm 

Da, 



(3.17) 



depends still on T = tj — U, i.e. M = M(T). As short-hand notation we shall use in the 
following 

<x '\U(T)\x>=^ J Daexp{£rf dt a (t) ■ V 1 ■ a (*)} < X K (T)\ X > , (3.18) 
M = J Daexp!^^J dt a (t) ■ v' 1 ■ a , 

and 

U a (T) = T[exp|-^y dth ff (t)}] 
with single-particle Hamiltonian 

h (a(t)) = h a {t) = (K + a(t)) ■ p . (3.19) 
The time-ordering operator T is necessary since in general 

[M*),M*')1 + °- ( 3 - 2 °) 

In (3.18) we have represented the time-evolution operator of a system with time-independent 
two-body interaction as superposition of time-evolution operators with time-dependent single- 
particle Hamiltonian h a (t) defined through some collective field <r(t). The superposition, being 
written as functional integral, contains the typical Gaussian weight factor in the field a. 

4 Ambiguities of the auxiliary field 

There are two sources of ambiguity or, positively put, of freedom in the choice of the auxiliary 
field a (t) [3, 9]: 

1. There are various ways to divide the Hamiltonian into one- and two-body parts, 

H = T + V = (T + U) + (V - U) = H + H' , (4.1) 

a fact well-known from the static shell model. All such versions are equivalent in an exact 
treatment, however, approximate schemes of calculating physical quantities will obviously lead 
to different results. 

With the choice (3.4) for the density operator p, the mean field a, determined in lowest order 
of the stationary phase approximation to the functional integral (3.18), is of Hartree type. The 
Fock term appears only after quadratic corrections are taken into account. In contrast, choosing 
Pa s = a a a $ as density operator in (3.5), one obtains in lowest order the Fock-term only which 
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in turn is corrected in second order by the Hartree-term. With = a) a and A<5 7 = a 7 as 
it is possible to introduce pairing fields in lowest order. 

All told, the solutions a of the stationary phase approximation just define a starting point for 
higher orders of the stationary phase approximation. Their convergence rate will, of course, 
depend on the choice of the starting point. 

2. The second ambiguity lies at the core of the functional integral representation introduced in 
section 3. In the expansion of the second term in the exponent of (3.12), 



{% 1 % e 2 

~h ed ' v ' p \ = 1 ~ ' v ' p ~ 2F ^ ' v ' p ^ '" " ' 

the linear term does not contribute to the integral as it is odd. The contribution of (4.2) to 
the integral in (3.12) stems from the quadratic term in (4.2) which actually is of order e, since 
with the Gaussian weight factor the dominant values of a are of order e~2. We may, therefore, 
arbitrarily modify the coefficient of the term linear in a in (4.2), without changing the value of 
the integral (3.12). If we choose this coefficient as 

WafiyS = Vaf3~f6 ~ Va/35-y , (4.3) 

then we will have included both Hartree and Fock terms already in lowest order stationary 
phase approximation (SPA). 

This linear term is important when one determines the stationary value a° of the integrand 
in (3.12), but does not contribute to the integral (3.12) or finally the functional integral in 
a corresponding to (3.18). Conversely the a-quadratic terms in (4.2) do not influence the 
stationary value of a, however, they must be taken into account when the functional integral 
over the fluctuations around a° is calculated. 

5 Stationary phase approximation 

To evaluate the integrals in (3.18), we shall use the stationary phase approximation (SPA). We 
illustrate the method for a real, one-dimensional integral 

oo 

7(0 = J dxexp(-f(x)/l) . (5.1) 

— oo 

For small I the integral is dominated by the stationary points of f(x). If there is only one 
solution / (xo) = 0, then expansion of f(x) upto second order in (x — xo) gives 

1(1) = exp(-/(s o )/0 j/ dsexp (- {x ~ pi f" (x ) ± • • •) j (5.2) 



exp (- f(x ) 1 1) 



I 2nl 

rS5' 
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assuming a minimum (f"(xo) > 0). The complex variant of (5.1), (5.2), extended to the multi- 
dimensional case, will be applied to the integral (3.18). The role of the small parameter / is 
then taken by %. Hence the use of (5.2) implies a semi-classical expansion. 

Rewriting (3.18) as 

<X\U{T)\ X >= ^ J £>(rexp{is[<7(f)]} (5.3) 

with 

S = - J dta(t) ■ v' 1 ■ a(t) - ik\n< x'|T[exp(- ^ J dth a (t))]\x > , (5.4) 



the saddle point condition for the functional S, 

5S 



5 a ai (t) 

reads explicitly, in matrix notation, 



for all a, 7 , (5.5) 



, <x'\U a o(t f ,t)pU t7 o(t,t i )\x> 



W <X\Uao(tf,t i )\ X > 1 ^ 

with U a o (t" , t') the time-development operator under the single-particle Hamiltonian h a o (t). 
Note that a will in general be complex, although the a-fields were originally introduced as 
real integration variables. If we take x 1 X as Slater determinants and generate x as 

\X >= U a o(t f ,t i )\x> , (5.7) 
then (5.6) simplifies to the TDHF mean field 

a o / t N = < x\U*° (U, t)pU a o (t, U)\x > _ v , 5g . 

<x\x> 

which refers to the initial time ti only. 

The value of the functional integral (5.3) with action S from (5.4) is in lowest order 

<X\U(T)\ X >= ex p{^/ dta°(t)-v- 1 -<j°(t)}<x\UAT)\x> ■ (5.9) 

The norm Af(T) is canceled by the quadratic correction term corresponding to (5.2), if we 
ignore the second term in S, eq. (5.4). This term carries an i ft-factor which makes it weakly 
oscillating compared to the first term. This term, dropped in (5.9), generates the Fock term if 
one starts from the original version of equations (3.4) and (3.5) with Hartree term only, and it 
cancels the self-interaction term in (3.6). It also reproduces correlations of the random phase 
approximation [3]. 

Taking x, x as Slater determinants, the equation of motion for the density operator is obtained 
from (3.18), (3.19) and (5.6) by defining Slater determinants 

!*(*)>= U a o(t, U)\ X > ; <*(*)! =<x\U*o(t f , t) . (5.10) 
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The time dependence of the respective single-particle functions ipi (t), ipi (t) is then governed 
by 

in — \i/>i(t) >= h a o(t)\^i{t) > ; -ih— < ^{t)\ =< ^(t)\h a o(t) (5.11) 

with time boundary conditions 

iPi(U) = x * ; 4>i{t f ) = x\ • (5.12) 

For the density operator 

= f m ><^wi (5.13) 

it", < i:(t)\i:(t) > 

one then immediately confirms eq. (1.7), 

ih^-p = [Ko (t),p\ . (5.14) 

Equation (5.14) has to be solved under time-boundary conditions, at t — and tf, rather 
than as initial value problem like (1.1), (1-2). Each transition x ~^ X w ^ have its own mean 
field, and this allows the method to describe exclusive reactions in contrast to TDHF which at 
best can be used for inclusive reactions. It is also worth noting that TDMF needs two sets of 
single-particle wave functions, ipi (t) and ipi(t), while TDHF is formulated in just one set which 
develops from the initial state in the mean field u. 



6 Static approximation 

On the exact many-body level the connection between the resolvent and the time-development 
operator is simply a Fourier transformation 

<X\(E + IK-HY 1 \ X >= dTexpj^(£ + < X '\ exp \x > (6-1) 

o 

for infinitesimal k > and with x, x as initial and final states. On the mean field level the 
connection between the resolvent and time-development operator is more intricate due to the 
non-linear nature of self-consistent mean field methods. 



The basic assumption to be made is that a (t) and hence h a (t) vary slowly in time so that we 
may approximately neglect their time dependence. In this static approximation we have from 
(3.18) and (3.19) 

< x'\ exp {-^Ht) \ X >= J Da exp { JL( ff • • a) t] < X '\ exp {-1 K t) \ X > 

(6.2) 

so that (6.1) reads: 

<X\(E + H)- 1 \x>= (6.3) 
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ft J U v y J AA(T) 

In iV-particle Hilbert space with 



N 

K = E (<7) , (6.4) 

m = 1 

we can use the inverse Fourier transformation on the single-particle level, 

f i 1 i °[ ex P x'k Wm T f 
exp<^ - -(/i m - i/c m )T^ = — / dw m — with T>0 (6.5) 

— oo 

for particle m, to obtain an expression for the many-body resolvent in terms of single-particle 
resolvents g rn = (w m - h rn + m m ) _1 



< x \(E + i« 



H r^>= -s Q N J dT jkJ D °J (n (66) 



A' 



X 



exp + m -E^+2 (7 ^ 1 ' °") T } < x I II ( Wm ~ h m + iK m ) l \x> 



) m = l 

with the constraint k = k — J2 K>m > to ensure the existence of the T-integral. 

m 

The product of single-particle operators in (6.6) acting on the Slater determinant x (to the 
right or on x to the left) gives 

( II 9m) A\xi X2---Xn>= A\ (gi Xi) (92 X2) ■ ■ ■ {9n Xn) > (6.7) 

\m = l / 

since the antisymmetrizer A commutes with the symmetrical product TT 9m- Equation (6.7) 

m 

suggests to introduce two sets of single-particle wave functions 

I > = 9m I Xm > 5 < ¥ m \ = < X m I 9m (6.8) 

or equivalently 

(W m - h)\<p m (<X, W m ) > = | Xm > ; < Vm (°", W m ) \ (lV m - h) = < Xm I ■ ( 6 -9) 

Note that on the single-particle level we may now drop the label m on the Hamiltonian h being 
the same for all (identical!) particles. In this representation 

< X \(E + in - H)-'\ X >= -i (A) J dT J D, J {U dwm ) ^exp {i S [<r, »,„, T]} 

(6.10) 

with 

S = {e + ik+)-(j-v~ 1 -(t-Y^ w m ) T - in E m < Xm I Vm > ■ (6.11) 

\ m / m 
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In the last term in (6.11) we can replace < x' m I ¥m > by < (f' m | Xm > using the bra- rather 
than the ket-equation in (6.8). In both versions we assume, with the same arguments as in 
section 2, that the overlap matrices < x m I ¥n > an d < ip m \ Xn > are chosen diagonal. Eqs. 
(6.5) and (6.8), (6.9) establish a Fourier transformation on the single-particle level between 
i\) m (T) and tp m (w m ): 



Wm{T)> = 




dw m exp |- -w m T^ \tp m {w m ) > ■ (6.12) 

— oo 

Equations (6.9) have, apart from the normalization factors Aj, A^, the formal structure of (2.7). 
However, energies w m were introduced as real variables in contrast to the complex rji of (2.7), 
and h in (6.9) is defined in terms of the a-fields rather than through ip m , (p' m . Only after we have 
solved all integrations of (6.10) in SPA, we will recover the inhomogeneous TIMF equations 
(2.7) together with the correct value of the resolvent matrix element (2.9). 




7 Path and energy integrals 

These integrals will be solved simultaneously in SPA, followed by the T-integration. While 
the energies w m and auxiliary fields (J Q, ey £1X6 introduced as real, independent variables, their 
stationary values will be complex and coupled through self-consistency. 



The stationarity equations with respect to a and w m read: 

6_S 
Ia~ 



rn N 

' To ■ v' 1 - ih ]T - Am ! s -f 1 = (7.1) 



^ Xm 


6 

5a 


lfi m > 


^ Xm 


ip m > 



m = l 

#^ = -T -ih V < X ' m 1 Vm > = . (7.2) 

dw n ^1 <Xm\ l fm> 

The above derivatives of ip m can be calculated with the help of equations (6.9): Multiplying 

$ S h 

(w m - h) — \(p m > ~ t— I ¥m > = (7.3) 
da da 

by < (p' m | gives, with h from (3.19), 

5 



Explicitly: 



< X m I 7— I > = < I P I <Pm > • (7-4) 

o a 



< Xm I 7 I > = < <P m | 4 a l I > = < <P m | « >< 7 I > • (7.5) 

d a ai 
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Derivatives of (6.9) with respect to w n give 



Snm I | Vm > + (w m - K) \ if m > \ = . (7.6) 



dw 

Multiplying by < (p' m | and using again (6.9) leads to 



d 

< X m I ~ I > = - S mn < <P m \ <Pm > • (7.7) 

o w n 

We insert relations (7.4) and (7.7) into (7.1) and (7.2) to find the explicit stationarity equations. 
The stationary value of a is given by 

*° • tr 1 = % £ <v4Mv?m> (7.8) 

T m <X m \ l Pm> ^ J 

or 

a" = - V <^\P\^> - V . (7.9) 

T m <Xm\ ( Pm> K ' 

Explicitly with labels of single-particle states 

ff a 7 = y E / | «/»««7 = Y ^ <v' |g > (7 ' 10) 

m Am I J m ^ Am I 

0,6 

using the completeness of single-particle states /3, S referring to the representation (3.5) of the 
Hamiltonian H . The stationary value is determined implicitly through 

- % -T (7.11) 
or, using (6.9) again, explicitly by 

< - £ m (7.12) 



< 




> 




<^m 


> 



i ft 






> 




I 7 ~~ 


< 


V?m 


> 





with 



< <p' m 






< V?m 


V?m 


> 



(7.13) 



After T-integration in section 8, E 1 rather than T 1 will appear in relation (7.12). Combining 
(7.10) and (7.11) identifies a° as Hartree-Fock field, 

N ' i i 

" 7 = ^ |^m> ' ( } 

m=l rm I rm ^ 

since according to section 4 the matrix v is assumed to be antisymmetrized. A slight modifica- 
tion of (4.3) also allows to eliminate the unphysical self interaction in (3.6) so that K reduces 
to the kinetic energy operator. 
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The normalization factors in (6.10), Af(T) and (i / 2n) N , are eliminated by the Gauss correc- 
tion to the lowest order SPA which approximates the integrals by the integrand taken at the 
stationary point (cr°, w°). The Gauss correction to the cr-integral, 



det 



' 5 2 S ' 
5a 5a' 



= (det(Tu- 1 ))' 



(7.15) 



cancels M{T) according to (3.14), if the a-dependence of the ln-term in (6.11) is neglected. 
This ln-term carries a factor h as compared to the leading term. Hence the above assumption is 
justified in semi-classical approximation. For the Gauss correction to the w-integration, each 
integral can be treated separately as the integrand factorizes. From (7.2) and (7.7), together 
with (6.9) and (7.12), it follows that the leading term in the second derivative of S with respect 
to w m is given by 

(7.16) 



d 2 S 



— t 2 

h 



and the total correction factor is 



2nih 



^-[J \S"{w m ) i 



f 2n 



T 



N 



(7.17) 



The correction to (7.16) involves fluctuations of the single-particle propagators which go beyond 
the mean field concept. 



To this order the matrix element (6.10) of the resolvent reads 

D(E + ik) =<x'\(E + in - H)' 1 ^ > 




exp 1 (e+Io* 



(7.18) 

1 -<t o -^£^ + ^t}< x '|0>, 



after replacing w° m by e° m according to (7.12). With e / \J2 ir = 1.084..., the normalization factor 
(z / 2 7r) is not exactly cancelled by the Gaussian correction. However, this is just a constant , 
real factor which does not harm the analytical properties of the resolvent. 



8 T— integration 

With regard to SPA, we extend the range of integration for T to — oo, introducing the step 
function in its integral representation, 

G(T) = — J ds 6RST „ for k" > . (8.1) 
2%i J s — i k 

— oo 

Then, with e / \[2~k w 1, we have 

oo oo 

D(E + in) = ^—j- J dT J ds expj^S[s, T]J (8.2) 
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with 

S[s, T\ = (e+ \a° ■ v' 1 ■ a° - £ e° m + a + i^j T 

-ihJ2 ln <x' m \ l Pm> -ihNhi + ih \n(-s + i k") . (8.3) 

Simultaneous SPA to both integrals gives as stationarity equations 

9S ^ ih ih , ^ 

a7 = r + 737^ = ° » --'" = ^ < 8 - 4 > 

and 

f| = £ + ■ v- 1 ■ <r° - £C + s + in = (8.5) 



or - s = E + \a° ■ V 1 ■ a - "£e° m 



where «" are now superfluous. Eq. (8.5) is exact though it looks as if in (8.3) only the term 
proportional to T had been taken into account. The contributions from the remaining terms 
cancel exactly by virtue of eqs. (6.9). 

For the Gauss correction to the lowest order SPA of the T- and s-integrals, one has to calculate 
the Hesse-matrix M of S with respect to T and s since S is not additive in the variables T and 
s. One finds 

u= {Uv) (8 - 6) 

with 

U - ° < 8 - 7 > 

in semi-classical approximation. Hence detM = — 1, and the Gauss correction factor for the 
T- and s-integrations is ^ — (2nhi) 2 = 2nh, which cancels the corresponding factor in (8.2). 
To this order of SPA we have 

D(E + i K ) = exp{ls[T 0?So ]} = 1 <x'|0> (8.8) 

= (^ + ^°-^ 1 -^°-Ec) <X>>, 

with = 0(<7°,iy^j). The same result is obtained following an alternative line of reasoning: In 
the above semi-classical approximation, with S linear in T and derivative d S / dT = E + 
\<j° ■ v^ 1 ■ a° — J2 £° m + s independent of T, the T-integral reduces to a ^-distribution, 



m 



2 ^/^e X p{l(£ + ^ 

(8.9) 
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With (ih/T) from the stationarity condition (8.4), the remaining s-integral can then be solved 
exactly, 

5 ( E + \ a° ■ v- 1 ■ a° - E £° m + s 

ds^ 2 ^(_ S )^< X '|0 >= (- 8o ) N - 1 <x\<l>> (8-10) 

s 

— oo 

with s from (8.5), in complete agreement with (8.8). 

9 Conclusions 

With cr° from stationarity condition (7.14), the single-particle Hamiltonian h(a°), eq. (3.19), 
has precisely the structure of h of section 2 with expectation value 



< <p' m 


(* + 






(t + u) 




< v' m 






> 



< ( Pm\ ( Pr 

From (7.8), (7.11) and (7.14) one finds 

1 _ v -l _ 1 < I °"° I V 9 ™ > _ 1 ^ V^m V^n I u I ^ ^ ^ 



2 2 



m 







\v\<f m 


¥n > 


< V'm 1 




><<Pn 


<Pn > 








i h 


= r) = 




- £ ° 


= jr , 

-L 



so that with (7.12), (8.4) and (8.5) 

E+-a-v -a-}^e m = E — — — — — = rj = w m - e m = — , 9.3 

2 - < | > T 

and we can identify the (complex) stationary value with r^ m of eq. (2.8). Hence eqs. (6.9) 
and (2.7) differ only by factors A;, A- and we have to identify 

Am. 



We make now use of (2.13) to obtain the final result for (8.8) 

D(E + in) =< x\(j) > = < <f)'\x > (9.5) 

which fully agrees with equation (2.9) for the stationary value of functional F in the TIMF- 
method. In summary we can state that we have derived the inhomogeneous single-particle equa- 
tions and the stationary value of the Green function of the TIMF-method from the TDMF in 
static, semi-classical approximation, with the Gauss-correction to the lowest order SPA taken 
into account. 



Under various conceptual and practical criteria, TIMF can be placed between TDMF and 
TDHF (see also table 1): Both TDMF and TIMF methods calculate a mean field for each 
transition x ~^ X ■ Hence the respective transition amplitudes and resolvent matrix elements 
refer to exclusive reactions, while TDHF is applicable to inclusive reactions only. TDMF is 
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conceptually superior to TIMF, the latter being a static approximation of the former method. 
Alternately put, TIMF is "local" in energy whereas TDMF is "non-local" in time. From the 
practical point of view, TIMF requires solving inhomogeneous, complex equations of Hartree- 
Fock type for given energy. In this respect it is comparable to TDHF which is "local" in time. 
The problem of TDMF lies in combining self-consistency with given boundary conditions in 
time. No practicable algorithm for this highly "non-local" problem seems to exist for use in ac- 
tual numerical calculations. The practical advantage of TIMF becomes even more pronounced 
when considering the S-matrix. TDMF then has to calculate three auxiliary fields rather than 
one [2, 3] while the corresponding T-matrix of TIMF [6, 7] is obtained again from inhomo- 
geneous, complex equations of Hartree-Fock type as for the above resolvent, with only slight 
generalization of the inhomogeneities. 
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Mean Field Methods for Reactions 



Basic quantity 



Time-Dependent 



< x' I exp{-i(tf-U)H} \ X > 



Time-Independent 



Fourier 



Transformation 



< X '\ (E-H)- 1 \ X > 



Mean Field Approximation 

restriction of \' 

TDHF <= 

\x'>=u h {t f -ti)\x> 

Equation of motion 

[h, p] = ihp 

P(U) 

Type of problem 



TDMF 



[h, p] = ihp 
to be solved for given 
P(U) , p(t f ) 



initial boundary 

condition problem in t: 
"local" "non-local" 

Representation of density p in basis 



orthogonal 



ip m (t) > , < ip m (t) 
biorthogonal 



Single-particle Hamiltonian h 

Hermitian non-Hermitian 

Specification of asymptotic channels 

no yes 

Reactions described 

inclusive exclusive 



static 



approximation 



TIMF 



[K p] = v(p X ' - Px) 



inhomogeneous problem 
parametric in E: 
"local" 



> , < <p' m 

biorthogonal 



non-Hermitian 



yes 



exclusive 



Table 1 



